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ABSTRACT 

We show that the Weibel-mediated coUisionless shocks are driven at non-relativistic propagation 
speed (0.1c < V < 0.45c) in unmagnetized electron-ion plasmas by performing two-dimensional 
particle-in-cell simulations. It is shown that the profiles of the number density and the mean ve- 
locity in the vicinity of the shock transition region, which are normalized by the respective upstream 
values, are almost independent of the upstream bulk velocity, i.e., the shock velocity. In particular, 
the width of the shock transition region is ~ 100 ion inertial length independent of the shock velocity. 
For these shocks the energy density of the magnetic field generated by the Weibel-type instability 
within the shock transition region reaches typically 1-2% of the upstream bulk kinetic energy density. 
This mechanism probably explains the robust formation of coUisionless shocks, for example, driven 
by young supernova remnants, with no assumption of external magnetic field in the universe. 
Subject headings: shock waves — plasmas — instabilities — magnetic fields 
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1. INTRODUCTION 

CoUisionless shocks are important phenomena in the 
universe; they heat and compress interstellar or inter- 
galactic plasmas and often accelerate charged particles 
or cosmic rays. The dissipation mechanism of the coUi- 
sionless shocks is highly complex and generally involved 
with instabilities in the coUisionless plasma. 

The Weibel instability is driven in anisotropic or 
counter-strea ming coUisionless plasmas (^W eibel 1 119591 : 
iFried I Il959f ) and generates strong magnetic fields. 
CoUisionless shocks inevitably generate an anisotropy 
within their transition region, where the upstream 
and downstream plasmas are mixed up, therefore, the 
Weibel instability can develop and generate strong 
magnetic fields there if the ba c kgrou nd magnetic field 
is weak. iMoiseev fc Sagdeevi ()1963D considered such 
a process in non- r elativ istic coUisionless shocks and 
iMedvedev &: Loebl (|1999f l discussed the magnetic field 
generation by the Weibel instability in relativistic coUi- 
sionless shocks of gamma-ray bursts. 

In this context, the Weibel instability in relativistic 
counter-streaming plasmas has been investigated ex- 
tensively by means of particle-in-cell simulations (e.g., 
Kazimura et al. I [19 98': Gruzinov '2001"; Haruki & Sakai ' 

Nishikawa et al. 2003; 
Recent 
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Frederiksen et al. I l2004l :~ ISpitkovskv I 120051) . 



research with simulations showed that the "Weibel- 
mediated" coUisionless shocks can form at rela tivistic 
propagati on speed in ele ctron-positron plasmas dKato I 
[2007: Chang et al. II2008I) and also in electron-ion plas- 
mas (jSpitkovskv I l2008l ) without background magnetic 
field; in these shocks the dissipation mechanism is effec- 
tively provided by the magnetic fields generated by the 
Weibel-type instability in the shock transition region. 
In the universe, these shocks can be driven associated 
with relativistic phenomena, such as gamma-ray bursts, 
jets from active galactic nuclei, or pulsar winds. 

There is also a possibility that these Weibel-mediated 
shocks are driven even in non-relativistic phenomena. 
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e.g., shocks of supernova remnants whose propagation 
speed is typically ~ 1, 000 — 10, 000 km s~^. However, in 
this case, the Weibel instability can be inefficient because 
the linear growth rate of the instability is proportional to 
the flow velocity V . Therefore, it is not clear yet whether 
this kind of shocks can exist in non-relativistic regime. In 
this paper, however, we show that the Weibel-mediated 
coUisionless shocks propagating at non-relativistic speed 
can form in unmagnetized electron-ion plasmas by per- 
forming numerical simulations. This also suggests a pos- 
sibility to proof the formation of these shocks in labora- 
tory plasmas with large-scale laser facilities. 

2. METHOD 

In order to investigate the coUisionless shocks in 
electron-ion plasmas without background magnetic 
fields, we performed numerical simulations. The simu- 
lation code is a relativistic and electromagnetic particle- 
in-cell code with two spatial and three velocity dimen- 
sions (2D3V) which was developed based on a sta ndard 
method described by iBirdsall fc Langdon I (|1991h . The 
basic equations are the Maxwell's equations and the (rel- 
ativistic) equation of motion of particles. The simulation 
plane is the x—y plane. The z axis is perpendicular to the 
plane. Since we consider coUisionless shocks in unmag- 
netized plasmas, the electromagnetic fields are initially 
set to zero over the entire simulation box. The boundary 
condition of the electromagnetic field is periodic to each 
direction. 

In the simulations a coUisionless shock is driven by 
means of the "injection method." There are two rigid 
walls at the left-hand side (smaller x) and right-hand side 
(larger x) of the simulation box and these walls reflect 
particles specularly. Initially, both electrons and ions 
are uniformly loaded in the region between the two walls 
at the bulk flow velocity of V in the -f-x-direction. At 
the early stage of the simulations, the particles located 
near the right wall were reflected by the wall and then 
interacted with the incoming particles, i.e., the upstream 
plasma, via some instabilities. Such instabilities induce 
strong electric or magnetic fields that provide dissipation 
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mechanisni of the shock and, eventually, a coUisionless 
shock is formed. Note that the frame of the simulation 
corresponds to the downstream rest frame. 

3. RESULTS 

We have carried out a series of simulations with up- 
stream bulk velocities oiV = 0.45c, 0.2c, and 0.1c, where 
c is the speed of light. The thermal velocity of electrons 
is typically taken as one-tenth of the bulk velocity and 
that of ions is determined so that the temperatures of 
electrons and ions are same. Because of the computa- 
tional power, we used a reduced ion mass of rrii — 20me. 
Here, we take uj~^ as the unit of time and the elec- 
tron skin depth A^ = ctOp^ as the unit of length, where 

cjpo = (47rneoe^/TOe)^" is the electron plasma frequency 
defined for the electron number density of the upstream 
plasma, Ueo- We also define the ion inertial length, 
Xi = {■mi/rne)^''^Xe. The units of electric and magnetic 
fields are both taken as E^, ^ B^, = c(47rneo'Tie)^/^. 

3.1. Results for V = 0.45c 

The simulation for V = 0.45c was performed on a 
grid of N^ X Ny = 4096 x 512 with - 27 particles 
per cell per species. The size of the simulation box is 
Lx X Ly = 2240 Ae x 280 Ag. (We choose the grid size so 
that the system is st able against the cold beam n onphys- 
ical instability, c.f. iBirdsall fc Langdon I (|1991h .) The 
number density of ions obtained from the simulation at 
ujpet = 2100 is shown in Fig.[TJ The color shows the num- 
ber density normalized by the upstream value, namely, 
neo . The upstream plasma flows from the left to the right 
and goes through the transition region which has a fila- 
mentary structure and then reaches the almost uniform 
downstream state (x/Ag > 1900). 




Fig. 1. — Number density of ions for the upstream bulk ve- 
locity of y = 0.45c at ujpct = 2100. The left-hand side and 
right-hand side are the upstream and downstream of the shock, 
respectively. Filamentary structure is visible in the shock transi- 
tion region, 1300 < x/Xc < 1900. 

Figure [5] shows the time evolution of the number den- 
sity of ions. The number density shown in color is av- 
eraged over the y-direction. The transition region, or 
"shock front" , which is visible as a steep increase in 
the number density, propagates upstream at an almost 
constant speed (~ —0.18c measured in the downstream 
frame) after ojpet ~ 500. Note that the particles that 
are reflected at the right wall at the early stage of the 
simulation, which are remains of the initial condition, 
eventually fade away at later times. 

Figure [3] shows profiles of quantities averaged over the 
y-direction at tOpct — 2100: (a) the ion number density, 
(b) the mean velocity in the x-direction, and (c) the en- 
ergy densities of electric and magnetic fields. Both num- 
ber density and mean velocity rapidly change to reach the 
downstream values through the transition region which 
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Fig. 2. — Time evolution of the number density of ions for 
V = 0.45c. The horizontal and vertical axes are the x-coordinate 
and the time, respectively. The color shows the number density av- 
eraged over the j;-direction. The shock transition region, or "shock 
front" , is visible as a steep increase in the number density. 

extends with the width of VF ~ 500 Ae. It is remarkable 
that a strong magnetic field is generated within the shock 
transition region. The energy density of the magnetic 
field reaches approximately 1% of the upstream bulk ki- 
netic energy density (measured in the downstream rest 
frame), Uke = neo{'mi + me)V^ /2. This strong magnetic 
field deflects and isotropize the particles coming from up- 
stream; this provides an effective dissipation mechanism 
for this coUisionless shock. 




Fig. 3. — Profiles of (a) the ion number density, (b) the mean 
velocity in the x-direction both normalized by the respective up- 
stream values, and (c) the energy densities of the magnetic (solid 
curve) and electric (dashed curve) fields both normalized by the 
particle bulk kinetic energy density of the upstream plasma. 

Figure |4] shows the current density in the x-direction, 
Jx, and the z-component of the magnetic field, Bz, at 
ojpct = 2100. It is evident that a number of current 
filaments exist within the transition region. These fila- 
ments consist of ions which are surrounded by the diffuse 
electrons similar to the case of the relativis tic counter- 
streaming plasmas (jFrederiksen et aini2004f ). The fila- 
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mentary structure seen in the number density in Fig. [T] 
reflects the existence of them. As seen in the figure, 
these filaments generate the strong magnetic field around 
themselves which is also observed in Fig. [3] (c). This 
structure would be formed through the Weibel-typc in- 
stabi lity like the relativistic shocks i n unmagnetized plas- 
mas (iKato 1120071 iSpitkovskv II200I . 



Ca) 



- V = C.45n. 





Fig. 4. — (a) The current density Jx and (b) the magnetic field 
Bz in the vicinity of the shock transition region at cUpct = 2100. 



3.2. Dependence on upstream bulk velocity and ion mass 

We also performed simulations with the bulk velocities 
of 0.2c and 0.1c. In both cases we confirmed that the 
structure of the shock is qualitatively identical to that 
oi V = 0.45c case. The ion number density and the 
magnetic field obtained for V ~ 0.1c are shown in Fig. [51 




Fig. 5. — (a) The ion number density and (b) the magnetic field 
_Bz obtained from the simulation for V = 0.1c. 



It should be noted that these Weibel-mediated shocks 
have similarities in the profiles of some quantities under 
changing the upstream bulk velocity V. Figure [H] shows 
the profiles as Fig. [3] at different upstream bulk veloc- 
ities: V = 0.45c, 0.2c and 0.1c. Here, the horizontal 
axes are the x-coordinate normalized by the ion iner- 
tial length including the relativistic effect, A* = F^/^A^, 
where F = [1 — {V/c)'^]~^^'^ is the Lorentz factor of the 
upstream bulk velocity. It is remarkable that the pro- 
files in (a) and (b) are almost identical independent of 
V. In particular, the widths of the shock transition re- 
gion are W ~ lOOA* in all cases. The energy densities 
of magnetic fields within the transition region also reach 
typically 1-2% of the upstream bulk kinetic energy den- 
sity at maximum independent of V. 

In order to see the effect of the mass ratio, we per- 
formed simulations for mi /me = 50 and 100 with V = 




Fig. 6. — Profiles of (a) the ion number densities, (b) the mean 
velocities in the x-direction Vx, and (c) the energy densities of 
magnetic fields Ub for V = 0.45c (red solid curves), 0.2c (blue 
short-dashed curves), and 0.1c (green long-dashed curves). The 
horizontal axes are the x-coordinate normalized by the ion inertial 
length (including the relativistic effect). A*. The normalizations of 
the profiles are same as in Fig. [3] 

0.45c. The normalized ion number densities for the dif- 
ferent three mass ratios are shown in Fig. [7] as in Fig. [S] 
(a). The width of the transition region is approximately 
W ~ lOOA* in all cases and this suggests the width may 
be the same order even for the real mass ratio. However, 
there is a slight difference in the profiles and this will be 
investigated in detail in the forthcoming paper. 




Fig. 7. — Profiles of the ion number densities for different ion-to- 
electron mass ratios: rrii/m^ = 20 (red solid curve), 50 (blue short- 
dashed curve), and 100 (green long-dashed curve). The upstream 
bulk velocities are V = 0.45c in all cases. The horizontal and 
vertical axes are same as in Fig. |6](a). 



4. DISCUSSION 

The similarities in the normalized profiles shown in 
Fig. [S]can be explained as follows. The structure of the 
transition region is essentially determined by the coa- 
lescence history of the ion current filaments. Here, we 
simply approximate the coalescence of two filaments as 
that of two straight currents. In this approximation, the 
current in each filament is proportional to the bulk ve- 
locity y, and a dimensional analysis of the equation of 
motion of the filaments shows the time scale of the coa- 
lescence of two filaments due to the Lorentz force, r, is 
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proportional to V^^. Thus, if we consider that the width 
in which the scale of the filamentary structure becomes 
twice is given by AW ~ Vt, it is independent of V. On 
the other hand, the downstream end of the transition 
region is determined by the saturation of the magnetic 
field of the ion current filaments; i.e., the typical gyro- 
radius of ions becomes comparable with the filament size, 
Tg^i ^ rf, like the c ase of electron-positron current fila- 
ments (jKato l[2005| ). This gives the condition rf ^ Xi, 
which is independent of V. Therefore, if we consider the 
initial scale of the ion current filaments is on the order of 
the electron skin depth Ae , the total width of the transi- 
tion region is also independent of V. 

Though the results described in the previous section 
were obtained for unmagnetized plasmas, these Weibel- 
mediated shocks can also be driven in magnetized plas- 
mas if the magnetic field is sufficiently weak. One of 
the criteria for neglecting the upstream background mag- 
netic field, Bq, may be given by rg^i/W ^ 1, where 
Vg^i is the typical gyro-radius of ions defined for Bq 
and W is the typical width of the shock transition re- 
gion. If we define the magnetization parameter a as 
(T = [BQ/87r]/[neo(TOi + rne)V^/2] and use the result 
obtained in the previous section W ^ 100 A*, the cri- 
terion can be rewritten a <^ 10^*. For shocks in su- 
pernova remnants (SNRs) of ^ 1000 years old, typical 
parameters are Tieo ~ 1 cm~^, V ^ 2000 km s^^, and 



Bq ^ 3/xG. In this case, the magnetization parameter be- 
comes a ~ 10~^ — 10"^ and therefore there is a possibility 
that the Weibel- mediated shocks are driven in SNRs. For 
younger S NRs with higher shoc k speeds of > 10, 000 km 
s^^ fe.g.. ' Revnolds et al. Il2008[ ). these shocks would be 
formed more easily. 

If the Weibel-mediated shocks are driven even at the 
velocity of ~ 1000 km s~-^, they can be investigated 
in laboratories because laser experiments are capable 
of generating coll i sionle s s plasma flow at t his velocity 
(iNishimura et al. 1 119811 : iRipin et al.1 Il990t ). For this 
plasma flow with Ueo ~ 10^" cm~^ (even at this den- 
sity, the plasma can be collisionless during the shock 
formation), the width of the shock transition region is 
estimated to be VF ^^ 2 mm. This is a reasonable value 
as the plasma size produced with the present-day large- 
scale laser facilities. To demonstrate the shock formation 
in laboratories is a challenging attempt but it would be 
particularly worthwhile as o ne of the key subjects for 
laboratory astrophysics (e.g.. iTakabe 112004 ). 
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